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ABSTRACT 

We present radial velocities for 2 045 stars in t he Sm all Magellanic Cloud (SMC), 
obtained from the 2dF survey by 'Evans e t al.| ( 2004[ ). The great majority of these 



stars are of OBA type, tracing the dynamics of the young stellar population. Dividing 
the sample into ad hoc 'bar' and 'wing' samples (north and south, respectively, of the 
line: 5 = — 77°50' + [4a]', where a is in minutes of time) we find that the velocities in 
the SMC bar show a gradient of 26.3± 1.6 km s~^ deg~^ at a position angle of 126±4°. 
The derived gradient in the bar is robust to the adopted line of demarcation between 
the two samples. The largest redshifts are found in the SMC wing, in which the velocity 
distribution appears distinct from that in the bar, most likely a consequence of the 
interaction between the Magellanic Clouds that is predicted to have occurred 0.2 Gyr 
ago. The mean velocity for all stars in the sample is +172.0 ± 0.2 km s~^(redshifted 
by ~20 km s~^ when compared to published results for older populations), with a 
velocity dispersion of 30 km s^^. 

Key words: galaxies: kinematics and dynamics - galaxies: Magellanic Clouds - stars: 
early-type - stars: fundamental parameters 



1 INTRODUCTION 

By studying the kinematic and chemical characteristics of 
the Milky Way and nearby galaxies, past interactions and 
the history of mass assembly can be explored. The Mag- 
ellanic Clouds are of particular interest in this context as 
their dynamical evolution and star-formation histories ap- 



pear to be closely entwined (e.g., Yoshizawa & Noguchi 
20031. For instance, the Magellanic Stream has generally 



been thought to arise from tidal disruption by the Milky 



Way ([Gardiner fc Noguchilp96l [Putman et al]p98l [Con 



nors et al.|2006 1, although new proper motion results (Kalli- 
vayalil et al.|20 0^ suggest that the Clouds may be on their 
first passage about the Milky Way, calling into question the 
origins of the Stream ( [Besla et al. 12007] ) . 

From observations of the neutral- hydrogen gas, |Math-| 
ewson & Ford ( 1984[ ) advanced the concept that the SMC 
is comprised of two components, with a 'Mini-Magellanic 
Cloud' being effectively torn off by a past interaction with 
the LMC, leaving the remnant of the SMC. These com- 
ponents had been identified previously by |McGee fc New-] 
|ton| ( |1981| ), who also noted additional features at larger 
and smaller velocities. Observations of the H I gas at bet- 



ter spatial resolution ( Staveley-Smith et al. 19971 have re- 



vealed a much more complex structure, comprised of hun- 
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dreds of expanding shells, presumably driven by stellar-wind 
outflows and supernovae, which led the authors to describe 
the SMC as 'distinctly frothy and filamentary'. Subsequent 
observations of three large supershells were reported by 
[Stanimirovic et ah] ( [1999] ), with the most dominant of these, 
304A, potentially accounting for the twin components seen 
previously in the main body of the SMC. 

To date, the most comprehensive survey of stellar kine- 
matics in the SMC is the study of 2046 red giants by [Hams] 
& Zaritsky (20065, who concluded that this old population is 



supported by its velocity dispersion (rather than rotation). 
Previous, more limited, studies have focussed on luminous 
supergiants ( Ardeberg fc Maurice|1977 1979 Maurice et al 



[19871, Cepheids ( [Mathewson et al.[[198 8), and on the older 



stellar populations, e.g., carbon stars ([Hatzidimitriou et al. 



itriou et al.[[1993| . 



1997 Kunkel et al. 20001 and red-clump stars (Hatzidim- 



The proximity and significantly sub-solar metallicity of 
the Clouds have also made them popular targets for spec- 
troscopic studies of stellar evolution. This was the primary 
motivation for our own spectroscopic survey of the young, 



blue population of the Small Magellanic Cloud ( Evans et al, 
[2004[ hereinafter Paper I), which resulted in spectral clas- 
sifications for 4161 stars, dominated by early (OBA) types. 
Here we present radial velocities for 2045 objects from that 
sample, covering a relatively wide spatial area (including 
the 'wing' region) and providing kinematic information for 
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a component of the SMC population which has previously 
gone largely unexplored. 



2 STELLAR RADIAL VELOCITIES 

Our spectroscopic survey used the 2-degree-field instru- 
ment at the Anglo- Australian Telescope (2dF at the AAT). 
2dF was a multi-fibre instrument which allowed up to 400 
intermediate-dispersion spectra to be recorded simultane- 
ously, using two spectrographs ( Lewis et al.|2002 1 . The main 
observing programme consisted of 18 fields observed in two 
runs: 1998 September and 1999 September (fields 1-12 and 
13-18, respectively). The survey area is shown in Fig. [l] 

Our spectra all cover at least the wavelength range 
~3900-4800A (with about half including H/3) at a resolving 
power of -R ~ 1600 {Av ~ 190 km s~^), with signal-to-noise 
ratios in the range 20-150 (see Paper I for details). While the 
data were not optimized for very precise studies of velocities 
for individual stars, the moderately large sample, and ex- 
tensive spatial coverage, encouraged us to investigate global 
trends in the overall kinematics. 

The 4161 observed targets cover a moderately wide 
range in spectral type (with classifications from mid 0-type 
to, in a few cases, late G-type), and hence also in temper- 
ature and in spectral morphology. Although we compiled 
radial-velocity measurements using semi-automated proce- 
dures, we found that, as a consequence of this range and of 
variations in data quality, 'manual' measurements gave bet- 
ter results. Compelled to use strong lines because of signal- 
to-noise constraints, we measured heliocentric velocities for 
H/3, H7, H(S, H8, and H9 (He was excluded due to blending 
with Ca 11 H). It is these results that we report here. 

Many of the spectra show emission in the lower mem- 
bers of the Balmer series, precluding reliable measurement 
of the line centres. In the majority of cases this emission is 
probably nebular in origin, but in some it will be intrinsic 
to the star. Affected spectra were excluded from the outset, 
leaving 3369 stars with three or more 'good' lines. 



from the two spectrographs differ by more than 30 km s ^, 
leaving results from 11 fields (3, 5, 6, 8, 9, 13-18) with a 
mean offset (Specl-Spec2) of +3.1 ± 14.1 (s.d.) km s'^ (A 
cut at 20 km s~^ further excludes fields 8 and 9, but has no 
material affect on any results.) 

We also used measurements of 66 stars from ancillary 
spectra obtained in 2002; although these were taken in only 



2.1 External sources of error 

Our data are free from one common source of radial-velocity 
error, namely mis-centring of the target in the slit, because 
of 'scrambling' of the signal by the fibre runs. However, the 
two 2dF spectrographs were mounted on the top-end ring 
of the AAT, bringing other external factors, such as fiexure, 
into the radial-velocity uncertainties. Such factors are well 
known to 2dF users (though not well documented), with a 
signature of obvious problems being offsets in the velocity 
zero-points between the two spectrographs (e.g., the unpub- 
lished technical report by Stanford & Cannon 200^. 

Velocities from the two spectrographs appear well 
matched in our data throughout 1999 (as judged from com- 
parison of the median velocities from the two sets of different 
targets observed simultaneously at a given telescope point- 
ing). However, the 1998 observations show poorer agree- 
ment, with offsets of up to ~100 km s~^ between spectro- 
graphs (Fig.[2]|. We rejected all fields where median velocities 
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one spectrograph, the observations (described by Harries 
et al.|2003 l were designed for radial-velocity measurements, 
and repeat observations show a scatter consistent with in- 
ternal errors alone. This gives a final sample of 2045 stars, 
for which results are tabulated in Appendix A. 

2.2 Internal uncertainties 

Data given in the Appendix include the standard error on 
the mean of the measurements for each star (i.e. cr/y'n, 
where n is the number of measured lines). The mean stan- 
dard error from 2045 stars is 11.1 km s~^, which provides 
an estimate of the purely internal uncertainty of our veloc- 
ities. Only ~3% of stars have standard errors greater than 
25 km s~^, and these can generally be attributed to low 
signal-to-noise spectra, or to unrecognized nebular contam- 
ination. 

Velocities are available from two separate 2dF observa- 
tions for only 24 targets from the 2045-star sample; results 
are shown in Fig.[3](which includes, for completeness, results 
for an additional 14 stars with observations from rejected 
fields). The pairs of observations are distributed among 14 
field/spectrograph combinations, so that it is impractical 
to try to establish separately any systematic offsets from 
these data. Instead, we simply compute the mean absolute 
difference between observations of each target; 12.4 ± 2.1 
(s.e.) km s~^. It can be shown that the implied 1-a uncer- 
tainty on a single observation is a factor 0.886 smaller; i.e., 
11.0 ± 1.9 km s~^ This 'internal error' estimate incorpo- 
rates any contributions from plate-to-plate zero-point shifts 
and astrophysical variability in radial velocity, in addition to 
stochastic uncertainties in the measurements; the excellent 
agreement with the value given in the preceding paragraph 
therefore implies that stochastic uncertainties dominate. 

2.3 External checks 

Of independent results at our disposal, the high-resolution 
{R ~20,000, Av ~ 15 km s~^) VLT-FLAMES sur vey of 
massive stars in two SMC fields by |Evans et al.|(2006[ ) offers 
both the largest overlap with the 2dF data, and the highest- 
quality data. A comparison between 2dF and FLAMES ve- 
locities for 28 apparently single stars is shown in Fig. |4] 
There is excellent overall agreement, with a mean offset 
of -7 km s"^(FLAMES-2dF). The standard errors of the 
FLAMES velocities are small (typically 2-3 km s~^); if the 
errors add quadratically, and if none of the stars has intrinsic 
radial- velocity variations, the dispersion in differential veloc- 
ities (14 km s~^) provides an upper limit to the statistical 
uncertainty of the 2dF results of 13.5 km s~^. 

We have provisional velocities for a further 13 stars from 
a follow-up FLAMES study of the SMC cluster NGC 346, 
at lower resolution {R ~7,000, At; ~ 42 km s~^; Evans et 
al., in preparation). Results for these stars are are in rea- 
sonable agreement with the 2dF values (Fig. [4|, excepting 
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Figure 1. The approximate boundaries of tiie 2dF SMC survey, superimposed on a y-band image by M. Bessel (left; contours show H I 
column densities), and on an H I velocity map (right; small circles show fields observed by [Harris fc Zaritsky [[2006^ . The dense optical 
region constitutes the SMC bar, with the Wing being the lower-density eastern region. Adapted, with permission, from |Stanimirovic| 
|et al.| | [2004l l. 



four 09.5-Bl stars which have FLAMES velocities that are 
rather smaller than 2dF values (by ~40-70 km s"'^); these 
stars lead to a relatively large dispersion in the velocity dif- 
ferences. It is not implausible that these could all be single- 
lined binaries (even though there is only a 1 in 8 probability 
of all four stars having velocity offsets in the same direction) . 

Differential velocities are also shown in Fig. |4] for 18 
A-type supergiants/bright giants observed with the AAT's 
RGO spectrograph at R ~2,700 ( [Evans fc Howarth||2003| ) ; 
in these relatively cooler, less massive stars, effects such as 
nebular contamination and binarity would not be expected 
to be as important as in the OB samples. The 2dF mea- 
surements show a ~20 km s~^ offset (in the same sense as 
found for the FLAMES comparison), which is only ^10% of 
the combined resolution elements of the two datasets. If we 
suppose that the dispersion in velocity differences between 
the 2dF and RGO spectra arises from these datasets in in- 
verse proportion to their resolutions, the implied statistical 
uncertainty on 2dF measurements is ~15 km s~^. 

We also found three other apparently good-quality ve- 
locities in the recent literature for seemingly single stars in 
common with the 2dF sample. These are AzV 170 ([Walborn 



et al. 



2000 



2004p 



2dFS#1086), AzV 215 ([Trundle et al 

2dFS#1352) and AzV 235 ( [Evans et al.|[2004[ 2dFS'#1416, 
where '2dFS' numbers are from the catalogue of Paper I). 
Results are included in Fig. |4] 

The results of all these velocity comparisons are summa- 
rized in Table [l] All the indications are that the l-cr statisti- 
cal uncertainties in our 2dF velocities are ~ 13 ± 2 km s~^. 
There is some evidence that the 2dF results may be system- 



2 [Trundle et~aL] ( [2004[ also [Trundle fc Lennon|[2005| label tab- 
ulated velocities as 't^isr'! hut values are actually heliocentric; 
Dr. C. Trundle, personal communication. 



Table 1. Comparison of 2dF radial velocities from repeat ex- 
posures, and with those from external sources (discussed in Sec- 
tion |23]l. Tabulated velocity differences are in the sense 'source' 
minus 2dF. All targets used in the comparisons are OB stars, 
except for the RGO-spectrograph sample which is comprised of 
A-type stars. 



Source 


R 


N 


A{Vr) 


cr 


Median 










[km s-i] 




2dF 


1600 


24 




13 




FLAMES (HR) 


20 000 


28 


-7 


14 


-9 


FLAMES (LR) 


7000 


13 


-25 


21 


-20 


RGO 


2700 


18 


-18 


17 


-20 


'Other' 


20 000 


3 


-19 


(4) 


-18 



atically too positive, by perhaps ~10 km s~ , but there is 
no suggestion that this offset, if real, is a function of time 
or position. Given the general level of internal consistency 
demonstrated by Figs. |3] and |4] we are confident that our 
results are, at the least, internally robust. 



3 RESULTS 

Figure [5] shows the spatial distribution of our sample by 
velocity bin. By inspection, there is a trend for more positive 
velocities to the south and east in our sample; in particular, 
the wing region appears to have systematically more positive 
velocities than the bar. Arbitrarily dividing the sample stars 
into 'bar' and 'wing' according to whether they lie north or 
south of a line: 



-77°50' + (4a)' 
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Figure 2. Measured radial velocities from survey fields 1, 6, 10 & 13. Solid circles are velocities from spectrograph 1, open circles those 
from spectrograph 2. Fields 6 and 13 show no important systematic offset; Field 1 and 10 show larger differences (median velocities differ 
by 47 and 81 km s"'^, respectively), and were excluded from the analysis. 



(Fig. [sj where a is the right ascension in minutes of time) , 
histograms of the 2dF velocities for the two regions are 
shown in Figure [5] Also shown are the results grouped by 
spectral type; B-type spectra are grouped as either 'Early 
B' (earlier than B5) or 'Late B' (B5, B8 or B9), with the 14 
peculiar AF-type stars included with the FG stars. 

The most striking feature in Figure [6] is the difference 
between the velocity distribution of the early B-type stars 
compared with the later types. The early B-type distri- 
bution is near-symmetrical and single-peaked, centred on 
175 km s~^. This is in strong contrast to the FG stars which 
appear to have a double-peaked distribution, with maxima 
at 150 and 180 km s~^. This internal comparison reveals dif- 
ferent velocity components in the different spectral groups, 
i.e. with agcrl the relatively small number of 0-type stars 



^ The appearance of the distributions is, of course, a factor of the 
adopted bin-size. For convenience in Figure|6]we chose 10 km s~^. 
Adopting a bin-size of 15 km s~^ still reveals a clear difference be- 
tween the early B- and FG-type results, although the two maxima 
are less pronounced in the latter distribution. 



are also peaked at ~175 km s~^. The late B-type results ap- 
pear to have a similar distribution to the FG group, whereas 
the picture is less clear for the A-type stars. 

At first glance the 2dF results appear at odds with the 



velocities of OB A stars from Maurice et al. ( 19871, in which 



a double-peaked structure is seen, with the high-velocity 
peak much larger. However, their velocities (which include 
the measurements by Ardeberg fc Maurice||1977 1979 1 are 
predominantly from late B- and A-type stars which, when 
compared to our results in Figure [6] could account for the 
double-peaked distribution. Moreover, the remainder of the 
[Maurice et al.| sample comprises 23 O- and early B-type 
stars, of which 17 have positive Local Standard of Rest ve- 
locities (consistent with our O- and early B-type distribu- 
tions in Figure |6|, giving the larger, high- velocity peak in 
their results. 

One of the more interesting structures reported by 
IStanimirovic et al.| ( [19991 12004[ ) is the supergiant sheU (SGS) 
304A. This spans most of the northern region of our obser- 
vations, with almost 800 of the 2dF stars with radial ve- 
locities within the ellipse traced out by the shell (Table 2, 
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Figure 3. Differences in radial velocities from repeat observa- 
tions with 2dF. Open symbols indicate stars from rejected fields 
(Section [2TTJ; the outlying filled symbol (2dFS 2102) may indi- 
cate a spectroscopic binary. The horizontal dotted line is drawn 
at 13.0 km s~^, the internal 1-a error estimated from these data 
(Section p^. 
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Figure 4. Comparison of radial velocities from observations with 
high-resolution VLT-FLAMES spectra (filled circles, [Evans et al.| 
|2006[ ); low-resolution VLT-FLAMES spectra (open circles, Evans 
et al., in preparation); RGO spectra (open squares, [Evans fc] 
[Howarth] p003 ) ; other published high-resolution spectra (filled 
squares; see Section|2.3[for sources). 



Stanimirovic et ar]|1999 l. Over half of the stars are O- and 
early B-type stars, and their velocity distribution is indistin- 
guishable from those in Figure |6] The velocity diagram for 
SGS 304A shows tvyo distinct components at a pproximately 
130 and 170 km s~^ ( Stanimirovic et al.|[l999 l, presumably 
the origin of the offset components reported by [McGee fc] 



Newton ( 1981 1 and Mathewson & Ford ( 1984 1. There is sig 



nificant Ha emission at the edge of SGS 304A (Figure 11, 
[Stanimirovic et al.|[T999l ), with N66 (NGC346, the largest 
H II region in the SMC) at the northern limit. Indeed, the 
higher- velocity component is dominant in the R. A. -velocity 
slice at --72°13'40" (Figure 6, [Stanimirovic et al.|[2004| . It 
appears that the majority of recent star formation in this 
region, as traced by the 2dF early-type stars, is associated 
with the higher- velocity component. 



The FG stars observed with 2dF must be intrinsically 
luminous for them if at the distance of the SMC. Detailed lu- 
minosity classification of these later- type spectra was beyond 
the scope of Paper I, however at this point it is of interest 
to consider their physical luminosities, in order to estimate 
the range of ages and initial masses probed by the 2dF re- 
sults. The median magnitude of the observed FG stars, the 
majority of which are classified as F-type (187 from 229), is 
B = 15.9. Adopting {B — V)o = 0.37 and a bolometric cor- 
rection of 0.11™ for a F5 supergiant (e.g. [Johnson| 1966]), to- 
gether with an SMC reddening of E{B - V) = 0.09 (|Massey 



et a l.|1995 l and a distance modulus of 18.91 (Hilditch et al. 
2005), gives Mboi ~ —3.5. Evolutionary tracks from IChar- 



bonnel et al. ( 1993 1 then suggest an initial mass of ~5Mq , 
(corresponding to a main-sequence spectral type of ~B2-3) 
and an age of ~90-100 Myr. These are only rough estimates, 
but they demonstrate that the FG stars are beginning to 
probe a population significantly older than that traced by 
the early-type, main-sequence spectra. 



3.1 Mean velocities 

The mean velocities and weighted dispersions, by spectral 
type, are summarized in Table [2] The OBA-star velocities 
are significantly more redshifted than results both for our 
relatively small 'late-type' (FG) sample, and for many other 
older-star samples in the literature, which cluster around 
~150 km s~^ (see the summary in Table 3 of Harris & Zarit- 
[sky"2006). The latter difference arises in part because our 



spatial coverage extends further to the east than many pre- 
vious studies, which, coupled with the velocity gradient we 
have found, introduc es a positive bias. In particular, [Har-[ 
ris & Zaritsky (20061 report v — 145.6 ± 0.6 km s~^ for a 



red-giant population with a mean position about a degree 
west of that for our sample (Fig. [l|, where the mean veloc- 
ity of the OB A stars is ~160 km s~^. While the residual 
difference is certainly within the combined systematic un- 
certainties in the two datasets (and is of the size and sign 
discussed in Section [2.3[ ), it is of interest to note that the 
mean of the 7 velocities determi ned f or 50 eclipsing-binary 



OB stars by [Harries et al.| ( |2003| p] and [Hilditch et al.| ( |2005 l 



is 4-196.2 ± 24.8 (s.d.) km s"\ Their sample is, roughly, 
spatially coincident with [Harris fc Zaritsky[ s, but is a full 
50 km s~^ more redshifted on average. 

The difference for the OB stars is in the opposite sense 
to that which might be expected if the Balmer-line veloci- 
ties were affected by stellar- wind outflows (which should, if 
anything, give rise to a blueshifted bias). It remains to be 
seen how much of this is a consequence of prosaic effects 
(e.g, unrecognized line blends, or systematic offsets in 2dF 
spectroscopy), as compared to truly different dynamics for 
the massive-star population. 



3.1.1 Velocities in the Kl region 



The Kl region lies to the east of the SMC bar (de Vau- 
[couleurs fc Freeman|[1972| ; it is included in our wing sub- 

* Heliocentric corrections were applied in the wrong sense by Har- 
ries et al. (2003; cf. Hilditch ct al. 2005); we have corrected their 
7 velocities by -|-11 km to allow for this. 
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Figure 5. Spatial distribution of the 2045 stars witti radial velocities from the 2dF survey. The upper-left panel shows the whole sample, 
with the black line dividing the sample into ad hoc subsets of wing and bar stars (Sectionjsji. The axis of the maximum velocity gradient 
of the bar stars, for a centre at Olh 00m —73° 00', is shown in green (light grey in the printed version). The other panels show the 
positions of stars in 50 km bins. 



Table 2. Weighted mean velocities and Icr standard deviations 
of the 2dF sample, grouped by spectral classification. 



Type 


N 


V 


a 






[km s"-*-] 


[km s~-^] 


O 


83 


174.1 ± 0.9 


41.2 


B 


1251 


175.2 ± 0.2 


32.1 


A 


482 


170.9 ±0.4 


33.2 


FG 


229 


160.8 ±0.5 


35.1 


All 


2045 


172.0 ±0.2 


33.6 


Bar 


1572 


167.4 ±0.2 


33.3 


Wing 


473 


189.5 ±0.4 


28.6 



3.2 Global gradients 

To provide a simple quantification of the global velocity gra- 
dients in our dataset (Figure [5| we took two approaches. 
First, as a reasonably nonparametric treatment, we set up 
a series of spherical-coordinate systems, each defined by a 
right ascension and declination of zero latitude & longitude 
(increasing north and east, respectively), with the equatorial 
great circle rotated by position angle 6 (measured E from N) 
from the line of right ascension passing through the central 
point. For each such system we determined the coefficient 
of correlatiorj^ between stellar longitude and radial velocity; 
where the correlation coefficient r reaches a maximum, the 
maximum velocity gradient lies along the equator. 

Denoting the position angle of maximum velocity gra- 
dient as Smax, we find that the null hypothesis of no velocity 
gradient can be formally rejected with extremely high con- 



sample. Th is region con tains the N83 and N84 emission g^^^^^ 99.999%) at 0^,^, for any chosen central point. 



complexes (Henize 19561, and has been noted in the past 
as having a larger velocity of recession than both the bar 
and wing populations (e.g. [Ardeberg fc Maurice|1979l |Mar- 
tin et al. 1989[ ). We also find larger velocities for the 2dF 
stars in this region: v = 199.4 ± 1.7, a = 35.6, compared 
to V = 189.0 ± 0.4, a = 28.1 for the remaining 446 stars 
in our wing sample, where we have adopted a radius of 7.5' 
centred on N84A (a = Ihl4m3 7s, S = -73°18'27" 
|et al.|1989[|Bica fc Dutra|2000| ). 



Martin 



with 

6'n,ax ^ 126° + 0.24(50 - a) 

(where a is again the right ascension of the adopted central 

^ We computed the Pearson product-moment correlation coeffi- 
cient r, Spearman's p statistic, and Kendall's r; all give sensibly 
the same results. 
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point in minutes of time, and the optical centre of the SMC 
is near a = 50m) - roughly perpendicular to the major axis 
of optical (and H l) emission. The uncertainty in Smax is of 
order 5° (see below), while the small variations in r(6max) do 
not usefully constrain the location of the dynamical centre. 

As a more straighforward alternative, we also fit veloc- 
ities as a simple linear function of position, 

V = vo + AAa + BA5, 

where Aa, AS are co-ordinate offsets in degrees, measured in 
cylindrical equidistant projection. This is just a plane, with 
?;o the velocity of the plane at the origin of the co-ordinate 
system. Using the full dataset gives a velocity gradient of 
14.8 ± 1.0 km s"^ deg-^ at PA = 127.3 ± 4.6°, with vo = 
-tl73.47 ± 0.75 km s^^ at I'^.O, -73°. 0. 

Figure [T] shows velocities in a coordinate system centred 
on Olh 00m, -73° 00' along PA = 126°. The velocity gradi- 
ent here averages 14.2 km s~^ deg"'^ globally, but evidently 
most of the velocity change is across the SMC bar. Adopting 
the previous definition of 'bar' and 'wing' samples, the 473 
stars in the wing show roughly constant velocity (averag- 
ing +189.5 km s~^), while the bar stars show a gradient of 
26.3 ± 1.6 km s"^ deg"^ at PA = 126.2 ± 3.9° about a mean 
velocity of 167.4 km s~^. Note that this result is robust to 
the details of the adopted bar/wing separation, as shown in 
Table H 

The amplitude of global velocity changes in our data, 
~50 km s~^, is comfortably larger than any plausible source 
of error we have identified. However, although the H I results 



reported by Staveley-Smith et al. ( 1997 1 and Stanimirovic 



Figure 6. Velocity histograms of the 2dF results for the bar and 
wing sub-samples, and as a function of spectral type. 



|et al.| ( |1999[ ) also show more positive velocities in the wing 
than the bar (i.e., to the SE), the neutral-gas velocity field 
displays an overall NE/SW dominant gradient. The appar- 
ent discrepancy between the NW/SE gradient found here 
and the SW/NE gradient in the neutral-gas map can be at- 
tributed in large part to the region of largest H I velocities 
lying outside our survey region, around Ih 20m, — 71°30' 
(Fig. [l]). This region has a low density of stars (and gas), 
but purely in terms of H I velocity overwhelms the trend of 
increasing velocity to the SE which dominates our data, and 
which is also present in the H I measures. 

If we examine our results along the line of maximum 
gradient identified by [Stanimirovic etaTlpOOil their Fig. 3) 
we reproduce their results qualitatively, finding a gradient of 
10 km s~^ deg~^ (along a great circle at PA = 60°, through 
a = Oh 45m, S — —73°; Fig. [7|, though this still appears to 
be somewhat smaller than the neutral-gas gradient. 



3.3 Induced gradients 

An extended source moving perpendicular to the line of sight 
at its centre is not moving perpendicular to the line of sight 
elsewhere; the line-of-sight projection of the transverse mo- 
tion is therefore non-zero other than at the centre. Expressed 
in terms of observable quantities, the resulting induced ve- 
locity gradient is 

8.273 X 10"^ km s"^ deg"^ 

mas yr^^ kpc 

for proper motion /j, and distance D. 

From observations with the Hubble Space Telescope 
('ggr;, p<ariivayalil et ahlpOOel l determined a precise proper 




Figure 7. Left: projection of radial velocities about Olh 00m, —73° 0' at PA = 126°; stars assigned to the SMC bar (Section [3| are 
shown in green (right-hand grouping, light grey in the printed version), with those in the wing shown in red (darker grey). The solid lines 
show the gradient in the bar (26.5 km deg~^), and the relatively flat gradient further east, in the wing, as discussed in Section pi[ 
Right: projection of radial velocities about Oh 45m, —73° 0' at PA = 60° , the direction of maximum gradient in H I velocities (Sectionl3|; 
the line has a slope of 10 km s~^ deg~^. 



Table 3. Velocity gradients derived from 2dF velocities in the SMC bar, as defined by stars north of the lines given in the first column 
(where a is the right ascension in minutes of time). Small differences arc found depending on the adopted bar/wing separation, but the 
gradient and position angle are robust. 
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motion for the SMC. Such is the interest in their results that 
Piatek et al.|(2008) have recently reanalysed the ifST data. 



In Table|4|we list the induced velocity gradients arising from 
these two measurements, and also for the proper motion re- 
sults of' Kroupa fc Bastian| ( |1997| l and |Irwin| (| 1999*). In cal- 
culating the induced gradients we adopt D = 60.6 ± 1 kpc 
HHilditch et al.|2005| ). 

The proper motions yield induced velocity gradients in 
the same general direction as the gradient found in the 2dF 
results, but their magnitudes are too small to explain our 
results. Astrophysical and induced gradients will add vecto- 
rially, so that the implied magnitude of the residual astro- 
physical component of the early-type stars is ca. 18 km s^^. 



Harris & Zaritsky ( 2006 1 found a velocity gradient for 
red giants of 8.3 km s~ deg~^ - about the same magni- 
tude as the induced gradient inferred here from the results 



of Kallivayalil et al. (20061, but in an almost orthogonal di- 



rection (PA = 23°; that is, roughly along the major axis of 
the bar). Independent data are required to establish if the 
absence of induced gradients in red-giant velocities is a re- 
sult of errors in the measurements (velocities and/or proper 
motions) or is astrophysical (implying rotation of the red- 



giant population about the major axis, in the opposite sense 
to that found here, to compensate geometrical effects). 

3.4 Velocity dispersions 

If we take out the large-scale velocity trends in the data by 
subtracting a simple quadratic fit from the results shown in 
Fig. [t] (left-hand panel), the weighted l-a dispersion of the 
full set of velocity measurements is reduced from 33.6 km s~^ 
to 30.4 km s~^. Weighted dispersions for 'wing' and 'bar' 
stars are 28.6 and 33.3 km s~^, respectively; subtracting a 
linear trend from the bar sample reduces the latter figure 
to 30.7 km s~^, which is not significantly different from the 
wing value (as judged by an F test). 

Adopting 30.4 km s~'^as the most representative fig- 
ure for the observed velocity dispersion, then, assuming that 
observational and astrophysical contributions sum quadrat- 
ically, the intrinsic astrophysical velocity dispersion is ca. 
~28 km s"\ Although Fig. [t] shows that the total veloc- 
ity range across the observed area of the SMC is substan- 
tially larger than this velocity dispersion, most probably that 
range reflects tidal effects, not rotation. Thus while our ob- 
servational result contrasts in this respect with [Harris fc] 
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Table 4. Induced velocity gradients from published proper motion results. 
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Figure 8. Probability distribution function of radial velocities 
after correction for systematic gradients (cf. Section [3. 4| . A gaus- 
sian characterized by cr = 30.4 km s~^ is shown for comparison. 



4 SUMMARY 

We have presented stellar radial velocities measured for 2045 
stars from our 2dF survey of the SMC. We find a local veloc- 
ity dispersion that is approximately gaussian, with a ~ 30, 
similar to recent results for red giants (cr ~ 27.5 km s~^; 
[Harris fc Zaritsky| |2006[ ) , but the global velocity range 
is significantly larger, with the SMC wing receding some 
~20 km s~^ faster than the bar. We find a continuous dis- 
tribution of velocities across the SMC, but those in the wing 
appear to have a distinct distribution compared to those 
in the bar (left-hand panel of Fig. [7| , presumably a conse- 
quence of the last close-interaction between the SMC and 
LMC ^ 0.2 Gyr ago (|Murai fc Fujimoto]|1980| |Yoshizawa"fc 



Noguchi|2003l [Connors et al.|2006[ ). Within the bar we find 
a mean gradient of 26.30 ± 1.6 km s"^ deg"^ at PA = 126° 
(cp. 8.3 km s~^ deg~^ at PA = 23° for the red giants), 
although, as the H I maps show, this simple parameteriza- 
tion of the data almost certainly disguises a more complex 
dynamical behaviour. 



[Zaritsky| s result for red giants, we concur that the system 
is primarily supported by the velocity dispersion. 

Interestingly, the velocity dispersion for the 0-type 
spectra (Table [2| is found to be larger than that for the BA 
stars (at greater than 99% confidence). Weak nebular and/or 
wind emission could contribute to this, but only one of the 
83 0-type stars with measured velocities is commented as 
having weak emission at H7 (cf. Table Al from Paper I). It 
is plausible that undetected emission could be affecting the 
results, however Balmer emission is as common in the B-type 
2dF spectra as in the O-type stars (Table 7 from Paper I), 
so this does not suggest an explanation for the different dis- 
persions. We also note that the great majority of the 0-type 
stars are either dwarfs or giants, in which wind effects in the 
higher-order Balmer lines will be relatively minor. 

We suggest that the larger dispersion likely originates 
from undetected massive binaries and runaway stars. The 
latter suggestion relates to the interesting issue of the rel- 
ative velocities of field versus cluster stars (aside from the 
separate issue that 4% of 0-type stars appear to form in iso- 
lation, [de Wit et al.[[2005[). Using t he catalogue of clusters 
and associations from Bica & Dutra ( 2000 1, 18 of our 0-type 



stars with measured velocities are in the field (as defined us- 
ing the 30 pc criterion from [Massey et al.|1995[ ). For the 18 
field stars the (weighted) mean is « = 179.3 ± 2.2 km s~^ 
(cr = 48.9), as compared to t; = 173.2±0.9 km s~^ (cr = 39.3) 
for the 65 others. Unfortunately, with such a limited sample, 
these results are not statistically significant, but hint at a 
potentially interesting result, reinforcing the case for a more 
optimised kinematic survey of the Clouds. 
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APPENDIX A: CATALOGUE OF STELLAR 
RADIAL VELOCITIES 

The full catalogue of radial velocities for the 2dF sample 
is available through the on-line edition of Monthly Notices, 
and at the Centre de Donnees astronomiques de Strasbourg 

m 



(CDSQ 



Table [AT] shows an extract from the catalogue, to illus- 
trate the format. For each entry we give coordinates, adopted 
B magnitudes (including their source) and spectral classifi- 
cations, all from Paper I. In the final three columns we give 
the mean and standard error (Columns 7 & 8, in km s~^) of 
the measured velocities, and the number of lines used (Col- 
umn 9). 



ftp://cdsarc.u-strasbg.fr/pub/cats/J/MNRAS/XXX/YYY 
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Table Al. An illustrative section of the radial-velocity catalogue. Stellar radial velocities (vr) and the standard error on the mean 
(a/y/n) are given in units of km s"'^; the final column reports the number of Balmer lines that were measured for each star. Sources 
of adopted B magnitudes are |Massey|p002} , OGLE ( [Udalski et al.|1998| l, [Zaritsky et al.| | |2002| l, and 2dF (i.e. values from the scanned 
Schmidt plates, see Paper I), coded M, O, Z, and 2 respectively. 
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